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SUMMARY 

An earlier investigation of flame phenomena at the University of Michigan 
correlated flame speed with a non-dimensional parameter, the Peclet number. 
This investigation was undertaken to study the effect of inert diluents on flame 
speed and to determine whether the above correlation is valid when diluents are 
introduced into the unburned gases. 

A Bunsen burner flame was employed in this investigation using propane 
and ethylene as fuels. In addition to normal air, four volumetric nitrogen- 
oxygen ratios were used, ranging from 1.0:1 through 5.75:1. Data was re- 
corded at pressures from one atmosphere down to 1/4 atmosphere, using burner 
nozzles from 1/8 to 3/4 inches in diameter. For a given fuel and nitrogen- 
oxygen combination, an endeavor was made to hold constant the following 
quantities: oxygen-fuel ratio, jet velocity, and flame shape. 

As a result of this study, it has been concluded that: 

1. The previous correlation of flame speed with Peclet Number holds for 
the range of diluent inert gases investigated. 

2. Flame speed varies inversely with the concentration of the diluent gas. 

3. Apparent flame speed diverges more rapidly from true flame speed, at 
reduced pressures, as the nitrogen concentration is reduced. 

4. The slope of the true flame speed envelope varies inversely with nitrogen 
concentration. 

5. By application of adiabatic flame temperatures, it is possible to predict 
flame speeds throughout the range of diluent concentrations and pressures, 
knowing only the flame speeds at any two diluent concentrations. 


This study was conducted at the University of Michigan, Ann Arbor, Michigan 
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SYMBOLS 
Area, ft. | 
Activation energy, BTU/mol °R, 
Nitrogen - oxygen ratio, by volume. 
Specific heat, BTU/Ib °F. 
Diameter of burner, ft. 
Thermal conductivity, BTU/sec. ft. °F. 
Pressure, lbs/ft* or centimeters of mercury. 
Peclet number, dimensionless. 
Volumetric flow rate, ft?/min. 
Universal gas constant. 
Adiabatic flame temperature. 
Adiabatic flame temperature, pure oxygen. 
Gas stream velocity, ft/sec. 
Apparent normal flame speed, ft/sec. 
True normal flame speed, ft/sec. 
Apparent flame speed, ft/sec., as computed from the adiabatic 
flame temperature approach, eq. (5). 


Density at standard temperature and pressure, lbs/ft®. 


Subscripts 
B At burner. 
std Standard temperature and pressure, 14.7 psi, 520°R. 





INTRODUCTION 

A previous study conducted at the University of Michigan? demonstrated 
that the flame propagation rate in propane - air and ethylene - air Bunsen 
flames could be correlated with a non-dimensional parameter, the Peclet 
number. The writer reviewed some of the basic theory of the flame propa- 
gation phenomenon and explained the mechanism by which burner size and re- 
duced pressure change the apparent flame speed from the true flame speed 
which would be observed for a very large burner at a given pressure. He 
demonstrated that the apparent flame speed for different sized burners converges, 
at high Peclet numbers, on an envelope of true flame speed over the range of 
pressures investigated. This envelope of true flame speed appears as a straight 
line when plotted on a graph of flame speed versus the logarithm of the pressure. 
The relationship between Peclet number and flame speed was expressed mathe- 
matically as: Ya, /Vg = 1- e~PeP where 6 is a constant depending upon the 
reactants. It was concluded that: (1) the normal flame speed for the fuels 
studied varied inversely with the logarithm of the pressure from atmospheric 
pressures to as low as one inch Hg. absolute, and (2) the Peclet number must 
be kept large to avoid considerable reductions in flame speeds at sub-atmospheric 
pressure. 

The purpose of the present investigation is: (1) to determine whether the 
foregoing conclusions remain valid when the diluent nitrogen concentration in 
the combustible mixture is varied, and (2) to make a quantitative survey of how 
the flame propagation rate in the previously mentioned hydrocarbons is affected 
by the diluent nitrogen. For this purpose five nitrogen - oxygen mixtures were 


burned with each of the fuels, the minimum diluent nitrogen being 50% and the 


maximum 85% by volume in the "air'' reacting with the fuel. 

The work was conducted at the University of Michigan in the Aeronautical 
Engineering propulsion laboratory, East Engineering Building, during January, 
February and March 1956. It was performed as part of the third year curricu- 
lum in Aeronautical Engineering of the U. S. Naval Postgraduate School , 
Monterey, California, as partial fulfillment of the requirements for a Master 
of Science degree. It was financed by the Bureau of Aeronautics, Navy De- 
partment, Washington, D. C. 

We would like to express our appreciation to Professor R. B. Morrison 
for his overall guidance and enthusiastic interest in the project. Professor 


T. C. Adamson, Mr. R. E. Cullen, and Mr. M. L. Gluckstein also con- 
tributed many helpful suggestions which are greatly appreciated. To the 
Department of Chemistry and Metallurgy at the University of Michigan 
go our thanks for their generous assistance in the field of mass 


spectroscopy. 





EQUIPMENT AND PROCEDURE 

The apparatus herein described was provided by the University of Michigan, 
having been employed in previous combustion studies. A line diagram, pre- 
sented in Fig. 1, illustrates the configuration as employed in this investigation. 

The burner assembly consisted of a stainless steel tube approximately two 
inches in diameter, 24 inches long, and one-quarter inch thick. The com- 
bustible mixtures were introduced at the bottom of the tube through a three- 
quarter inch flexible line. The top of the tube was internally threaded to pro- 
vide for interchangeable bummer nozzles with nominal inside diameters from 
one-eighth inch to 1 1/4 inches. Screens just below the burner nozzles insured 
a constant low-level of turbulence in the flow and assisted in preventing flash- 
backs into the fuel lines. The upper four-inch portion of the burner tube was 
water-jacketed to provide the burner with as near a constant temperature as 
possible. Ignition of the mixture was accomplished by a high voltage, movable 
electrode at the burner nozzle. 

The pressure chamber, within which the burner assembly was mounted, 
was one foot in diameter and three feet high. A 100-centimeter mercury ma- 
nometer measured the pressure within the chamber. Sub-atmospheric pressures 
were maintained at desired levels with a Nash-Hytor vacuum pump, Manual 
control over the pressure level was exercised through a series of bleed lines 
and control valves in the vacuum lines. Two 1" x 7" diameter plate glass 
windows permitted visual inspection and photographing of the flame. 

A box camera was mounted on the pressure chamber at one of the windows. 
After obtaining optimun conditions of focus and parallax, the camera was 


locked into position to restrain further motion. Royal Pan 4" x 5'' sheet film 





was employed, with exposures ranging from f4 to f6.3 at 1/25 to 1/50 of a 
second. 

Because of the explosive nature of premixed combustibles, the fuel and 
nitrogen-oxygen supplies were stored in separate tanks. The flows were 
metered through individual rotameters, calibrated at standard conditions, and 
were united immediately downstream. Complete mixing was obtained by flow 
of the gases through several inches of tightly packed steel wool, this device 
also acting as a flame trap in case of flash-back. 

For any given fuel and oxygen-nitrogen mixture, data was recorded at 
pressures of 1, 3/4, 1/2, and 1/4 atmospheres. For each run, fuel-oxygen 
ratio and nozzle velocity were held constant while attempting to maintain 
similar geometric flame shapes. The fuel-oxygen ratio was that which gave 
maximum flame speed and was determined experimentally at atmospheric 
pressure prior to the run for each mixture. Although the fuel-air ratio for 
maximum flame speed varies slightly with pressure, this effect is secondary, 
according to Ref. 2. The nozzle velocity was that which provided a nearly 
constant geometric flame shape throughout the range of pressures and burner 
Sizes. At each pressure level, the flame was photographed and the following 
data recorded: rotameter flow rates, pressures and temperatures; chamber 
pressure; and room temperature and pressure. According to Ref. 1, the temper- 
ature of the unburned gases at the burner lip varies approximately 10 degrees, 
with negligible effect upon the results. In the calculations, this value was taken 
as constant at 85° F. 

The photographs were enlarged by factors of 10 to 40, depending upon 


original flame size. The average normal flame speed was computed from the 





equation 
Ve. = — (1) 


where Qp is the volumetric flow rate of the unburned gases at the burner lip 
and A, is the area of the inner reaction zone of the flame. This area was computed 
from the Theorem of Pappas: 

Ap = 277Lx (2) 
where L is the length of the inner reaction zone from base to tip of the photo- 
graphed flame, and x is the centroid of L with respect to the longitudinal axis 
of the flame. 

This investigation was limited to two fuels, ethylene and propane, the 
former being an unsaturated double-bond hydrocarbon and the latter a saturated 
single-bond hydrocarbon. The volumetric ratios of nitrogen to oxygen used with 
each fuel were 1.0:1, 1.8:1, 3.76: 1, 4.74: 1, and5.75: 1. These mixtures 
were prepared by the method of partial pressures and verified with mass- 
spectroscopic analyses. The ethylene, propane, nitrogen and oxygen used in 
this study had a guaranteed purity of 99% or greater. 

A serious limitation in the execution of this study was the inability of the 
rotameters to measure accurately very low rates of flow. In order to maintain 
constant burner velocity, flow rates decrease with decreasing pressure. With 
the simaller ateers. the desired flow rate was too low for the smallest rotameter 
available, and it was therefore impossible to obtain data at the lower pressures 
with these burners. Attempts to rectify this condition by employing larger 
burners failed in many cases because of rotameter limitations in the high ranges 


at the higher pressures or because of the onset of turbulence in the flames. 





This limitation might have been prevented had facilities been available for the 
use of premixed fuels, thus eliminating the necessity of two separate rotameter 


systems, both required to meter a wide range of flows. 





RESULTS AND DISCUSSION 

Due to the limitation described under EQUIPMENT AND PROCEDURE, the 
results in many regions were inadequate for presentation in their original form. 
In addition, considerable effort was expended in obtaining data at the lowest 
pressures. Although these points are of interest, they are not as critical in 
the results as are those at high Peclet numbers, i.e., high pressures. A 
multitude of points in these latter regions could have been obtained, with less 
emphasis on those at low pressures, had this effect been realized at an early 
date. As a consequence, it was impossible to arrive at a simple plot of flame 
speed versus pressure, and obtain the envelope of true flame speed, as shown 
in Figs. 5 and 6 of Ref. 1. Since this envelope is of paramount interest in 
obtaining the relationship of flame speed to Peclet number, it was necessary 
to derive this envelope by two rather indirect processes. These procedures 
will be outlined in detail later in this discussion. 

The Peclet number is the dimensionless parameter which is defined below 
and is equal to the product of the Reynold's number and the Prandtl number: 
- Pp Via PB ony 
a kK 


Since Cp/K is relatively constant, it is felt that if the Reynold's number were 


Pe 


used instead of the Peclet Number, similar results would be obtained. 

A series of photographs have been presented in Figs. 2, 3, 4 and 5 to 
illustrate the effects of pressure, burner size, nitrogen concentration, and 
fuel upon Bunsen-type flames. The effect of pressure is shown in Fig. 2. In 
this series of photographs, which were enlarged equally, the fuel is ethylene, 


using the 1. 8:1 nitrogen-oxygen mixture. Burner jet velocity and fuel-air ratio 





were held constant. As pressure decreases, the flame becomes more diffuse, 
the primary reaction zone increases in thickness, and the luminous secondary 
reaction zone disappears (not shown). In addition, the flame begins a gradual 
departure from the burner base, culminated by blow-off at some reduced pressure. 
Flame size gradually increases, indicating the retardation of flame speed, as 
calculated by eq. (1). 

The effect of burner size is shown in Fig. 3. The fuel, nitrogen concen- 
tration, jet velocity, and fuel-air ratio are constant, identical with those of the 
previous figure. The comparison is between the 1/8" and 3/16" burners at 
pressures from 1/2 to one atmosphere. The primary effect here is the ability 
of the larger burner to maintain a nearly constant geometric shape and size as 
pressure decreases. The attendant characteristic of a greater tenacity for the 
burner base may also be noted. These effects are attributed to the reduced effect 
of heat sink in the larger burner. The fact that the overall heat sink effect is 
reduced can be shown by considering that the heat sink is proportional to the 
diameter of the burner while the flame area is proportional to the diameter 
squared. 

In Fig. 4 is shown the effect of nitrogen concentration. Since burner size, 
jet velocity, and fuel-air ratio were variant it is only possible to surmise the 
effect of nitrogen from the photographs. However, burner size has a negligible 
effect here since no great changes have occurred in flame shapes between the 
3/4 and 1/2 atmosphere points. Likewise, although fuel-air ratios have changed, 
the fuel-oxygen ratios are within five percent. The greatest change occurred in 
jet velocity, decreasing from 13 to 2 feet per second as the nitrogen concen- 


trations increased from 64 to 83 percent. This indicates a large decrease in 
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flame speed, with decrease in reaction rate and flame temperature. The flame 
became more diffuse, with a reduced secondary reaction zone, and reduced 
tenacity. These effects are caused by the excess heat lost in raising the temper- 
ature of the added diluent from 85°F. to flame temperature. 

The effect of fuel is shown in Fig. 5. Here again jet velocity and fuel-air 
ratio are not identical for the two fuels. The fact that the ethylene reaction 
rate has increased over that of propane is indicated in the smaller flame sizes, 
with resulting increases in flame speed. This effect is attributed in part to the 
_lower activation energy of the ethylene fuel. It is also noted that the ethylene 
flame has greater tenacity and a greater secondary reaction zone (not shown). 

In correlating the results of this study with that of Ref. 1, it was desirable 
to obtain a relationship between flame speed and Peclet number; hence, a plot 
of apparent flame speed versus Peclet number was made for each of the fuels, 
using Cullen's data to establish basic trends, i.e., curves of constant burner 
size and curves of constant pressure. These curves are shown in Figs. 6 and 
7. The pressure curves are seen to reach a constant flame speed as Peclet 
number increases. This ultimate value of apparent flame speed is the true flame 
speed for that particular pressure and is the value needed to establish the true 
flame speed envelope on the flame speed versus pressure plots. 

From the trends thus established it was possible to draw curves through 
the few points in each of the Ng/Op runs and to make a fair estimation of the true 
flame speed for each point (Figs. 8 to 13). It is interesting to note that as the 
amount of diluent is decreased, both the curves of pressure and burner size 
steepen considerably and appear to peak at higher Peclet numbers. This results 


in a steeper slope for the true flame speed envelope as nitrogen content decreases. 
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The plots of apparent flame speed versus pressure, with the estimated true 
flame speed envelopes, are shown in Figs. 14 - 17. 

To take these results a step further, Fig. 18 is a plot of the ratio of View V¢ 
versus Pe. Although Ref. 1 combined the two fuels on the same graph, it is 
evident from these graphs of the several concentrations that different curves 
are obtained for propane and ethylene at each diluent concentration. Increasing 
diluent concentration results in a translation of the curve to the left. 

A second approach to this same problem, which gave similar results, 
involved consideration of adiabatic flame temperatures and activation energies. 
Due to the complexity of this approach, only the details involving propane as a 
fuel will be discussed. 

By the method outlined in Ref. 3, the adiabatic flame temperatures were 
computed for propane at each of the five nitrogen-oxygen concentrations. These 
temperatures were plotted versus nitrogen-oxygen ratio, as shown in Fig. 19. 
If it is assumed that this relationship is linear, then the equation relating the 
variables is 

Tg = Tg, - ke (3) 
where TE, is the flame temperature with pure oxygen, k is the slope of the curve, 
and C is the nitrogen-oxygen ratio. 

In Ref. 4, it was shown that 


Ae 





2RT; 
a 


where A, isthe activation energy, R is the universal gas constant and Ve 5 is the 


apparent flame speed. In this case, a single rate controlling reaction is assumed. 
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Substituting the relationship of eq. (3) into eq. (4) and letting a = A,/2R, 





a 
Tre. - kC 
O 
Miao Ee 
Then 
1 a 
hy ee 
ete TE, ~ kC Tg ~ kC' (9) 
Vy’ 0 
fo 


where Vf. is the apparent flame speed at any other nitrogen~-oxygen ratio, C’. 





Then, 
Ve 
oa 
V fa T_ ~kC Te -kC 
( Tz,-kC ) (Tg -kC') Ve 
etc ean Ma Ly pa (6) 
k(C' - C) vie 


By selecting values at a given Peciet number and burner diameter from the 
data of this investigation for Vfg, Vfa', C, and C', a value of "a" may be com- 
puted, which by definition, is constant for a given fuel, regardless of the value 
of C. Then by eq. (5), the apparent flame speeds for any nitrogen-oxygen ratio 
may be related to those for air, simply by application of appropriate values of C. 

Applying this procedure to the four non-air propane flames and plotting the 
results against logarithmic pressure produced the curves of Fig. 20. The curve 
for normal air is also included. Next, the envelope for true flame speed of the 
propane-air mixtures was taken from Fig. 5 of Ref. 1 and plotted. Using this 
curve in combination with those oe each nitrogen-oxygen ratio, it was possible 


to construct the curves of Fig. 21. 
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Comparison of the curves of Fig. 18 with those of Fig. 21 reveals a marked 
similarity. Since the means by which they were constructed were totally unre- 
lated, it is felt that the trends indicated therein are representative of the results 
to be expected from nitrogen dilution. The effect, as stated earlier is as follows; 
for increased nitrogen, the apparent flame speed approaches the true flame Speed 
at Peclet numbers lower than those for air, while for decreased nitrogen, the 
apparent flame speed approaches the true flame speed at Peclet numbers higher 
than those for air. 

Although it has become apparent that there exists a separate curve for each 
nitrogen concentration and each fuel, this fact does not invalidate the previous 
correlation of Peclet number with flame Speed, as proposed in Ref. 1. Each of 
the curves of Figs. 18 and 21 is still a close approximation to the general 


expression: 


SPor 
Vea/Ve - l-e 


as developed in Ref. 1. The variation is due to the “constant, & , which takes 
on new values for each fuel and each nitrogen concentration. Knowing the flame 
speed of a given fuel with air and with a second nitrogen-oxygen mixture, it is 
possible to obtain the value of f for this fuel and any other nitrogen-oxygen 
mixture, using the adiabatic flame temperature approach (eq. 5). 

For a number of reasons, it must be emphasized that these results are, 
for the most part, qualitative only. There were several sources of error in 
this investigation and the degree of error is unknown in some instances. As 
an example, consider the case of computing flame areas, where the most 
Critical area, the base of the flame, is the most difficult to determine due to 


loss of luminosity in the diffusion zone. Even with care and good judgment, 
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the determination is difficult, and probably best obtained by averaging the 
results of a number of separate runs. For the time allotted, this procedure 
was not acceptable, and it is therefore possible that an error of 20% could 
exist in some of the results presented herein. Other sources of error existed 
in the analyses of the nitrogen-oxygen mixtures, the flow rates indicated by 
the rotameters, particularly at the lower limits, and the assumption that the 
gas temperature at the burner remained constant at 85°F. Compared with the 
error incurred in computation of flame areas, these are all very small and 
therefore insignificant. 

A tabulation of the results is presented in Appendix J, and a complete 


set of calculations for a typical run is presented in Appendix I. 
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CONCLUSIONS 

For the range of pressures and nitrogen-oxygen ratios investigated, it is 
concluded that: 

1. The extensions of the classic Mallard-Le Chatelier model, as proposed 
in Ref. 1, is applicable to flames with diluent gases in ratios other than those 
found in air. | 

2. The apparent flame speed for the fuels investigated varies inversely 
with the amount of diluent present. 

3. The apparent flame speed diverges more rapidly from the true flame 
speed, at reduced pressures, as the nitrogen concentration is reduced. 

4. The slope of the true flame speed envelope varies inversely with 
nitrogen concentration; i.e., the true flame speed with pure oxygen decreases 
more rapidly than it would with air as pressure increases. 

5. By application of the adiabatic flame temperature approach, it is 
possible to predict flame speeds throughout the range of diluent concentrations 
and pressures, knowing aie the flame speeds at two different diluent concen- 


trations. 
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RECOMMENDATIONS 

Since the foregoing conclusions were made from observed trends rather 
than from a number of experimental points, it is felt that further experiments 
in this subject should be undertaken to provide more conclusive evidence. [If 
this should be done, the following recommendations may be of assistance: 

1. Use premixed fuel and nitrogen-oxygen concentrations to remove the 
necessity of metering separate flows. This will allow a wider range of data. 

2. For a given diluent concentration, obtain data at frequent intervals 
in the regions of high Peclet numbers for each of four or more burner sizes. 
This will determine the true flame speed envelope more precisely. 

3. Continue to record data at pressures of 1/4, 1/2 and 3/4 atmospheres 


to insure accurate plots in the regions of low Peclet number. 
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APPENDIX U0 
SAMPLE CALCULATIONS: 
The data in Table I were obtained in the following manner. Run #5 is used 


for the example. 


rotameter rdg. gage pressure Temp. 
No/Oo 20.90 9.6 psi 74°F, 
Propane 8.85 9.6 74 


Flow at standard conditions was obtained from rotameter calibration curves. 
= ie 

No/Oo9 0.1123 ft/min. 

Propane 0.00515 ft*/min. 

These flows were converted to the actual flow under standard conditions by 
a factor converting the pressure and temperature to standard, said factor from 
a graph in the rotameter handbook. 

Std. N2/O2 flow 0.1123 x 1.26 = 0.1418 ft?/min. 


Std. Propane flow 0.00515x 1.26 = 0.0065 ft/min. 


0.1483 fe int = Q std. 
't Air''-fuel ratio: .1418/.0065 - 21.8/1 


Oxygen - fuel ratio: Air 21% Oy by Volume 


21.8 x 0.21 


- 4.58/1 
; / 


Oo/F = 


Stoichiometric equation: CgHg t 502 —= 3CO, + 4H5O 
(O2/F) Stoich. = 5/1 
The Oo/F ratio of 4.58/1 used in run #5, the mixture which produced the 
best flame speed, was slightly richer than stoichiometric. 


From the A/F ratio the standard density of the mixture is obtained: 





Propane std. density 0.11616 ib/ft? 
Air std. density 0.0766 1b/ft® 


Mixture std. density: (0.11616 x 1)+(.0766 x 21.8) 
= 01saulb 


22.8 
From gas law: P = pRT 
14.7 x 144 
R = P#T = —————_ = _ 82.0 
.0783 x 520 


Assuming an average burner temperature of 545°R., the density of the mixture 


at the burner before combustion may be calculated; Pp = 14.38 psi. 








14.38 x 144 
f,= Pp/RTR = ————— = .0732 b/tt* 
52 x 545 
To get the actual flow at the burner, QB 
fstd Q, .1483 x .0783 
: | | «1585 ft? min. 
OB Qstd . 0732 
Velocity of gas at the burner, Up: 
Area of burner = 0.1188 in@ 
Ue ee eee «8. 2 «ft/sec. 
Ap .1188 x 60 
Apparent flame speed, Ven: 
Area of flame from photograph = 0.336 in? 
QB .1585 x 144 
Vea See ee C= s«12 13 ft/sec. 
Ag .336 x 60 
Peclet number: 
AR Via Dp Cp 


K 





Specific heat: 
O 
Cy ye BTU/Ib°F 


“ o) 
Cy propane = .576 BTU/Ib°F 


(A/F) by vol = 21.8/1 
21.8 x 28.9 
(A/F) by wt. = = 14.45/1 
1 x 44 


(14.45 x .25)+(1 x .576) 
ane $$$ $_$______—_—_—_—_—_———— = 0.262 BTU/Ib°F 
“ 15. 45 


Thermal Conductivity, K: 


K air - 424.5 x 107° BTU/sec ft°F 
K propane = 294 x 10° BTU/sec ft°F 


(14.45 x 424.5)+(1 x 294) 
K mi 10” «C= «416 x 10° BTU/SReeOR 
15. 45 


.O732 x 1.13 x .0324 x .262 


- 169 
416 x 1078 
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